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Abstract. A new procedure is proposed to generate planar 
asymmetric coronas with multi-streamer structures, as shown 
in SOHO/LASCO observations. The initial coronal magnetic 
field is separated into potential part and non-potential part. 
We first fit the potential part by using the sum of magnetic 
multi-poles, whereas the non-potential part is approximated 
by the magnetic field induced by some properly fitted current 
densities. Then the total field is numerically modified in terms 
of Maxwell’s equations. Based on the initial magnetic fields 
obtained in such a way, we can further compute various 
complex asymmetric corona structures by solving MHD 
equations. In order to verify this procedure, we compute the 
2-D coronal structures prior to the December, 1996 CME and 
August, 1999 CME events. The numerical result is on the 
whole in agreement with observations, and hence set proper 
backgrounds for further studying the propagation of CME in 
various coronal structures using numerical computations. 


1. Introduction 


The study of the coronal structure is an important topic in 
space physics and solar physics. It has been shown by 
observations that in most cases the corona has a variety of 
asymmetric multi-streamer structures. Since CMEs must pass 
the low coronal environment before moving into the 
interplanetary space, the study of coronal structure is 
fundamentally important, especially for quantitatively 
understanding the changes in space weather [Dryer, 1998; 
Feng and Wei, 1999]. There have been a great number of 
numerical works devoted to the study of the steady structure 
of the corona and CMEs. Pneuman and Kopp [1971] studied 
the 2-D steady structure of the corona by iteration method. 
Steinolfson et al [1982] simulated the steady structure of the 
corona by numerically solving the 2-D MHD initial boundary 
problem through relaxation method. Based on a symmetric 
coronal streamer background, Wu and Guo et al. [1997] 
investigated the formation and evolution of CME observed by 
LASCO (Large Angle Spectroscopic Coronagraph), and 
considered the interaction between flux ropes and a single 
streamer. So much attention is paid to the symmetric coronal 
structures, asymmetric multi-current sheet configurations, 
nevertheless, have frequently been observed in realistic 
coronal structures. Linker and Mikic et al. [1999] developed a 
3-D MHD model of the solar corona applicable to the Whole 
Sun Month time period, using measurements of the 
photospheric magnetic field as boundary conditions for their 
computation, and obtained results in good overall agreement 
with coronal and interplanetary structures. On the other hand, 
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since direct and accurate measurements of the coronal 
magnetic field cannot be made at present, some assumptions 
have been customarily introduced regarding its properties to 
make up for the lack of direct information. The measured 
photospheric field has been used as input to a potential 
field-source surface model to compute the coronal or 
interplanetary magnetic fields [see Linker et al., 1999, for 
example]. Zhao and Hoecksema have developed many 
methods [Zhao et al , 1993,1994,1999], such as the horizontal 
current-resource surface (HCSS) model, the horizontal 
current-current sheet (HCCS) model and the current 
sheet-source surface (CCSS) model, by respectively including 
the horizontal current, current sheet and volume current into 
the generating function © of the magnetic field or by 
adopting a good synoptic frame of the photospheric field 
corrected using the unseen surface adjustment to improve the 
potential field-source surface (PFSS) model. 

Due to the lack of measurement of coronal magnetic field, 
it is a very interesting problem how to construct the coronal 
structure (including the magnetic field and plasma properties) 
based on the coronagraph observation. This paper presents a 
new method to numerically simulate complex asymmetric 
corona with multi-streamer structures on the basis of the 
coronagraph observation from SOHO/LASCO. 

To verify the method, we apply it to reproduce the 
asymmetric multi-current sheet coronal structures prior to the 
December, 1996 CME event and August 1999 CME event. 
The structures obtained by the method are rather stable, and 
hence can be used as the background to further simulate the 
propagation and evolution of the CME mentioned above in 
the corona and the interplanetary space. 

The paper is organized as follows. We describe the 
scenario to generate asymmetric coronal structures with 
multi-streamer structures in section 2. In section 3, two 
examples are given to illustrate the method; in section 4, a 
brief discussion is made about the method. 


2. Description of the method 


In order to construct asymmetric corona with 
multi-streamer structures on the basis of the coronal 
observations, we proposed a magnetic field 
fitting-modification method. It consists of the following steps: 
1) fitting the initial magnetic field, 2) determining the other 
initial flow parameters, and 3) solving the MHD initial and 
boundary value problem based on the above initial field and 
flow parameters (using the nonreflecting projective 
characteristic boundary conditions). Step 1 is the core of the 
method, which account for its name. 


2.1 Fitting and modification of the initial magnetic field. 


When using the time-relaxation method to get realistic 
steady coronal structure, it is crucial to choose suitable initial 
magnetic field on the basis of the SOHO/LASCO 
observations. For this reason we develop the following 
method to fit and modify magnetic field. 
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It is well known that the magnetic field can be thought as 
consisting of two components: the potential part B, and the 


non-potential part B,. they satisfies 


Ve =0 
S (1) 
VxB,=0 
and 
VeB, =0 
= s (2) 
Vx B, = hJ 
separately. 


When the corona is in a simple and quiet state, its magnetic 
field can be taken as a potential field, otherwise, when the 
plasma interacting with the magnetic field, and the effect of 
the charge current cannot be ignored, the coronal magnetic 
field is not potential field at all. Both the potential and 
non-potential case are taken into account. In the following 
treatment of the magnetic field, we draw ideas from the 
theory of incompressible flow in hydrodynamics, and treat the 
magnetic field as treating velocity in the theory of 
incompressible flow, taking advantage of the analogy between 
magnetic field and the velocity in incompressible flow in 
hydrodynamics. 

2.1.A. Fitting of B,. Since B, satisfies equation (1), 
it can be represented as 

Bi =V¢ (3) 
then it follows 

Aġ=0 (4) 
that is, ø satisfies the Laplace’s equation. According to 
partial differential equation theory, in 3-dimension case, ¢ 
can be expanded as the following sum: 


(0) 
g (r) = 9(o)4 2 — 4.0 =) 
a) a XA r (5) 


where the right-hand side are successively dipole, 


quadrupole ...from the third term. 

In 2-D case, a similar expansion holds, if replacing 1/r by 
lnr. 

This suggests that we can fit the magnetic field B, by a 
delicate combination of multi-poles. That is, we place some 
dipoles, quadrupoles, etc, in appropriate positions, and 
compute the magnetic field strength of every such multi-poles, 
then take the sum of them as B, . 

2.1.B. Fitting of Bs Since B, satisfies (2), it can be 
approxiamted by the following way. 

1) Place some currents density 7,(x),J,(x), ,J,(x) in 


appropriate positions( in the current case, they should be 
perpendicular to the considered plane), and take their sum 


k 
SG) =>) (6) 

rl 

2) Use the formula 
ioo- te dTa (7) 
4n r 
to get the induced field, and take it as the approximation of 
Bp. 


By properly choosing the strength and position of every 
multi-pole and current density, we can construct an initial 
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magnetic field for a given corona. Starting with the initial 
magnetic field, we wish to further get the observed 
asymmetric coronal magnetic structure by time-relaxation 
method. 


2.2. Modification of B = B, + B, 


Generally speaking, the result obtained through 
time-evolution method starting from the above assembled 
initial magnetic field is not very satisfactory. It needs to be 
modified in the following way so as to get a better result. 

The following relations can be derived from equations (1) and 


(2) 


ð ôB 
a he) T + Ur (8) 
10°(rB,) af @ ay 
LO 22784) |= zru + Myr? Or (9) 
10°(rB.) O( @ a 
—~—— + =| r— (rB, ) |= -ur — 10 
r 00 +2(rZ08,)] Hrag (10) 
ff B,ds = 0 (11) 
{ Bydl = py! (12) 


here g ,B,, have the usual meanings, and / represent the 
current enclosed inside 1 Rs. Equations (11) and (12) represent 
the surface integration and the line integration along the inner 
boundary respectively. The reader should keep in mind the 
fact that we only consider the 2-D planar flow. 

We can use the above five equations to numerically modify 
the above fitted magnetic field as follows: 

Step 1: According to observed pictures, make proper 
modification on the component 8B, on the inner boundary so 


as to better approximate the observation. Usually, this can be 
successfully done by repeatedly comparing the computational 
result with the observation. Of course, educated experience is 
also important. Meanwhile, it should be noted that the 
modified B, has to satisfy equation (12). 

Step 2: Discretize equation (9) into difference form, use the 
B, values obtained in step 1 as the inner boundary condition, 
and take the normal derivatives of B, to be zero on the outer 
boundary. Then solve the difference equations by using the 
Gauss-Seidel iterative method to get the values of B, in the 
whole annual region from 1 Rs to 10Rs. 

Step 3: Determine the values of B, on the inner boundary 
by using equations (8) and (11). 

Step 4: Proceed as in step 2 to determine the values of 
B in the whole computational region from 1Rs to 10Rs by 
equation (10). 


2.3. Determination of other initial parameters 


Up to now we have described how to construct the initial 
magnetic field by using our fitting-modification method. 
Following the same way as in the simulation of axisymmetric 
coronal streamer, we can determine the initial distribution of 
other flow field parameters such as the density, velocity 
[Steinolfson. et al , 1982; Zhang and Wei, 1993]. That is, by 
solving the 1-D steady flow when the magnetic field is zero, 
we can obtain the distribution of radial velocity, density, and 
pressure along the radial direction, and take them as the initial 
values of the corresponding varieties. The latitudinal velocity 
is taken to be zero. 
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Figure 1. a) Observational picture of December. 1996 CME 
event. b) Computed density result at t=100.000s represented by b) 


brightness Figure 3. a) Observational picture of August, 1999 CME event. 
b) Computed density result at t=50.000s represented by 
. , brightness 
2.4. Numerical solution of the MHD equations 

Considering the interaction between the magnetic field 
and the coronal plasma, we apply MacCormack’s difference the above magnetic field as initial field and the other initial 
scheme to solve the time-dependent MHD equations, using parameters. The numerical solution after a sufficiently long 
time computation is taken as our steady coronal structure. The 
non-reflection projective characteristic boundary conditions is 
used as the boundary conditions on the inner boundary and 
outer boundary [Li and We1,2001]. 


3. Simulation of December, 1996 and August, 
1999 asymmetric coronal multi-streamer 
structures 


In order to verify the above procedure, we simulate the 
asymmetric corona prior to the December, 1996 and 
August,1999 CME events. The equations we use to describe 
the motion of the solar wind are the 2D planar ideal MHD 
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Figure 2. a) Initial magnetic field and velocity for calculating 

the corona prior to the December event: b) Magnetic and velocity Figure 4. a) Initial magnetic field and velocity for calculating 
at t=50,000s; c) Magnetic and velocity at t=100,000s; d) the corona prior to the August event.. b} Velocity and magnetic 
Magnetic and velocity at t=110,000s. field at t=50.000s 
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equations taking the gravitation of the sun into account. The 
computational domain is the entire annualar region from | Rs 
to 10Rs. The computational grid is 75 (radial direction) x 
288 (meridional direction). The steps along the radial 
direction increase as a geometric series, but the steps along 
the meridional direction are even. The initial temperature on 


the inner boundary is1.8 x 10° K, the maximal magnetic field 
is about 1.5 Gauss, y is taken as 1.1. 

Fig.1.a) shows the December, 1996 asymmetric coronal 
structure observed by SOHO/LASCO. It can be seen that 
there is a current sheet near the equatorial plane in the east 
side of the sun, and there are two current sheets in the west 
side, with one above the equatorial plane and another below. 
The whole structure is asymmetric. 

Fig.1.b) shows the result of density represented by 
brightness. The different values of brightness describe the 
changes of p — p,, where p is the simulation result of density, 


Pp is the initial density. The brighter part is where the larger 


difference of the density lies. 

Fig.2.a) shows the initial magnetic field obtained by our 
fitting-modification method and the superposed initial 
velocity distribution; Fig.2.b) to d) show the evolution of the 
magnetic field and the velocity. It can be seen that the flow 
basically reaches stable at t=100,000s and the velocity in open 
field regions is larger than that in closed field regions. The 
magnetic field result is on the whole in agreement with the 
inferred observation obtained by SOHO/LASCO. 

Fig.3.a) shows the August, 1999 asymmetric coronal 
multi-streamer structure observed by SOHO; b) shows the 
result of density represented by brightness(at t=50,000 
seconds). The location of the streamers approximately agrees 
with the observed one. 

Fig.4.a) shows the initial magnetic field for calculating the 
corona prior to the August event. 

Fig.4.b) shows the velocity and magnetic field at t=50,000s. 
Again we can see the basic consistence of magnetic 
configuration with the SOHO/LASCO observation. 


4. Discussion and Conclusions 


In the previous section, assuming that the December, 1996 
CME and August, 1999 CME events started from the 
asymmetric coronal structures as shown in Figure l.a and 
Figure 3.a respectively, we study both coronal backgrounds 
by using the procedure described in section 2. The 
computational results demonstrate the power of the 
fitting-modification method to simulate complicated coronal 
structures. More importantly, the method provides a way to 
numerically get rather complex 2-D coronal magnetic 
structures. Based on such coronal structures, we can 
quantitatively investigate the propagation of CME in a more 
realistic corona-interplanetary space. It should also be noted 
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that the method should be generalized to 3-D cases in order to 
get better results consistent with observations. Since the 
procedures described in section 2.1.A and 2.1.B obviously can 
be applied to 3-D case without any change, we know that 
there is no essential difficulty in such a generalization. 
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